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Abstract: We have investigated the head-on collision of a two-kink and a two-antikink 
pair that arises as a generalization of the model. We have evolved numerically the 
Klein-Gordon equation with a new spectral algorithm whose accuracy and convergence 
were attested by the numerical tests. As a general result, the two-kink pair is annihilated 
radiating away most of the scalar field. It is possible the production of oscillons-like con¬ 
figurations after the collision that bounce and coalesce to form a small amplitude oscillon 
at the origin. The new feature is the formation of a sequence of quasi-stationary structures 
that we have identified as lump-like solutions of non-topological nature. The amount of 
time these structures survives depends on the fine-tuning of the impact velocity. 
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1 Introduction 


Topological defects are of great interest due their ubiquity in many branches of physics 
like fluid mechanics, condensed matter, nuclear and particle physics and cosmology. The 
simplest form of topological defects are the one-dimensional kinks modeled by a scalar field 
(j) satisfying the Klein-Gordon equation, 


dV{4>) 

dt'^ dx‘^ d4> 


( 1 . 1 ) 


where V{(p) is the scalar field potential. A large variety of kink solutions arises for potentials 
with at least two distinct minima. The most common models of kinks [1, 2] are described 
by the potential V{(p) = (1 — /2 that defines the model, and V{(p) = 1 — cos((/>) 

known as the sine-Gordon model. 

The collision of kinks has been studied in details motivated mainly by the applications 
many areas of physics. The nonlinear nature of the problem governed by a time-dependent 
partial differential equation is responsible for new and unexpected features. We mention 
the fractal structure in the collision of a kink and an anti-kink in the cp'^ model [3-8], and 
also found in other models such as the parametrically modihed sine-Gordon [9, 10] and the 
(p^ models [11-13]. Remarkably, the exception is provided by the collision of sine-Gordon 
kinks in which the Klein-Gordon equation is solved exactly. 

Some years ago Bazeia et al [14] proposed a new class of topological defects in systems 
described by real scalar fields in (D, 1) spacetime dimensions. By restricting to D = 1 they 
have considered the family model with potential, 

v{cp) = ^^\r^/p-(pypf, ( 1 . 2 ) 

where the parameter p is related to the way the scalar field self-interacts. This model can 
be understood as a generalization of the <p^ model, which is recovered for p = 1. For p 
even, the case p = 2 is special and describes an unstable lumplike configuration [14]. For 
p = 4,6,.. the potential (1.2) requires that (/> > 0 or (/> < 0 under the change (p —>■ —cp. 
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As pointed out by Refs. [14, 15] it is possible to construct topological defect in the form 
(/>(■*■) = tanh^(x/p) {x > 0) and = — tanhP(x/p) {x < 0). Moreover, non-topological 
lumplike defects can be envisaged according to Ref. [16]. It can be shown that all defects 
with p even, topological or non-topological, are unstable and therefore of few interest. 

The new static structures arising when p = 3, 5.. connect the minima (j) = =tl passing 
through the symmetric minimum at 0 = 0 are called two-kinks defects [14]. As shown in 
Fig. 1, the two-kinks defects seem to be composed of two standard kinks symmetrically 
separated by a distance proportional to p. The energy density prohle reinforces this view. 
However, Uchiyama [17] presented for the hrst time a two-kinks structure after extending 
a hadron model. For the sake of completeness, the boosted two-kinks defects are described 
by the following exact solution. 




± tanh^ 


X — xq — ut\ 

J ’ 


(1.3) 


where ± denotes the two-kink [K) and the two-antikink {K), respectively. The parameter u 
stands for the velocity of the kink along the axis and xq the center of the two-kink/anti-kink. 
The total conserved energy E associated to the scalar held is calculated from. 


E = 



p^{x,t)dx, 


where the energy density is given by. 


(1.4) 




1 / 

2\dtj 




+ Vi^). 


(1.5) 


It can be shown that E = 4p/(4p^ — l)\/l — u'^ is the conserved total energy of the boosted 
two-kink / anti-kink. 

It is worth pointing out that the two-kink defects can describe the magnetic domain 
walls in constrained geometries [18, 19]. We also mention possible applications of two-kinks 
structures in the brane-world scenario [20, 21[, and in connection with processes of their 
formation in perturbed sine-Gordon models [22]. 



Figure 1. Profiles for the kink and its energy density for p = 3 and u = 0.4. The resulting 
structures are composed of two standard kinks symmetrically separated by a distance proportional 
to p [14]. 












The main purpose of this paper is to study numerically the collision of a two-kink-two- 
antikink pair with a code based on spectral methods. The paper is organized as follows. 
In Section 2 we have introduced the numerical method that can be applied to any one¬ 
dimensional scalar field model. We have performed numerical tests showing the accuracy 
and convergence of the algorithm considering the exact breather solution of the sine-Gordon 
model, and the collision of standard kinks of the cj)'^ model. In Section 3 we have proceeded 
by studying the head-on collision of a two-kinks pair for p odd. The outcomes are quite 
distinct from the collision of kink-antikink pair of the (j)'^ model. In Section 5, we have 
summarized the results and discussed future perspectives of the present work. 

2 The numerical method 

The Klein-Gordon equation (1.1) can be solved numerically using straightforwardly the 
Collocation method [23] as we describe in the sequence. The first step is to establish the 
approximate scalar field as, 

N 

(j)a{t,x) = '^ ak{t)il^k{x), (2.1) 

k=0 

where ak{t), k = 0,1,..,A^ are the unknown coefficients, the functions 'ipkix) constitute 
the basis functions chosen suitably as a generalization of a Fourier expansion, and N is 
the order of the series truncation. It remains, therefore, to compute the iV -|- 1 unknown 
coefficients afc(t) to complete the approximate solution. 

The solutions describing the dynamics of kinks are defined in the unbounded domain 
(—oo, -|-oo) and obtained after solving a nonlinear partial differential equation. In this vein, 
we mention two possible strategies for treating numerically this problem. The hrst is to 
consider the hnite difference scheme using staggered leapfrog integration [24] of the equation 
of motion in a hnite spatial domain with periodic boundary conditions. Pseudospectral 
method was used in a hnite domain with periodic boundary conditions for the (/>® model 
[12]. The second strategy which we will adopt here is to map the inhnite interval to a 
hnite domain. According to Boyd [23], it is possible to generate new basis functions for 
the inhnite interval as images under the change of coordinate of, for instance, Ghebyshev 
polynomials or Fourier series. There is a large variety of maps, some of them are more 
popular than others. In order to decide which map is more appropriate for the problem 
under consideration, we recall that a typical kink exact solution decays as powers of the 
hyperbolic tangent as \x\ —>■ oo. In this case, we have found that the logarithmic map 
[23, 25, 26] is well suited for the problem under consideration. The logarithmic map is 
given by. 


rj = tanh ^ ^ ^ ^ ^ (2.2) 

where r/ G [— 1,1], Lq is the map parameter and xq dehnes the origin of the computational 
variable rj. We can dehne the basis functions as the Ghebyshev polynomials, Tk{r]), or their 
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images under the map (2.2) defined on the physical domain x G (— 00 , 00 ): 


i)k{x) = Tk ^tanh Lq^^ )) ’ 

In what follows the collocation method is applied straightforwardly. The first step is to 
substitute the approximate scalar field into the Klein-Gordon equation to form the residual 
equation, 


N N 

ResKG{t,x) = ^ Vk{t)'ipk{x) - ^ ak{t)'4)l{x) 
k=0 k=0 



(2.4) 


where we have introduced the ‘velocities’ Vk = dk- Next, the unknown coefficients are 
determined with the condition of vanishing the residual equation at the collocation or grid 
points on the physical domain Xj, 


ResKG{t,Xj) = 0 , 


(2.5) 


for all j = 0,1, ..,N. The collocation points Xj are determined from the map (2.2) taking 
the collocation points defined on the computational domain, 

rjj = cos , j = 0,l,..,N. (2.6) 

Therefore, we ended up with a set of + 1 hrst order ordinary differential equations for 
the velocities Vk, which must be complemented with the set of + 1 differential equations, 
dj = Vj, to complete the whole set of equations for the unknown coefficients. While the 
first two terms of the residual equation can be expressed spectrally, that is, in terms of the 
unknown coefficients, the last term is treated using the physical representation. It means 
that instead of the spectral representation of the scalar field through the coefficients ak{t), 
we have considered the physical representation using the values of the scalar field at the 
collocation points, = (pait^xj). Both representations are related since, 

N 

= (paixj,t) = ^ ak{t)ijjk{xj). (2.7) 

k=0 

The resulting dynamical equations are, 


dk = Vk ( 2 . 8 ) 

Vk = Fk{aj,4’j), (2-9) 

where the first set relates the definition of the velocities and the second set arises after 
solving the relations represented by (2.5). The above dynamical system is complemented 
with A^ + 1 algebraic relations, 

4>j = 4>j{ak). ( 2 . 10 ) 
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To integrate the above equations one needs the initial data: 


</>(0,a:) =/(x), = gix), (2.11) 

from which the initial coefficients 0 ^( 0 ) and Xfc(O) can be calculated. 



Figure 2. Exponential decay of SEi-^s evaluated from t = 0 to tf = 70, after setting uj = 0.5 and 
to = 1.0 into the exact solution (2.12), and Lq = 10 for the map parameter. The saturation due to 
round-off error is achieved for N > 50. 



Figure 3. Exponential decay of SE,-^s for the collision of a kink and an antikink in model for 
u = 0.3. 


The classical sine-Gordon model offers exact solutions describing the interaction of 
kinks initially at rest or moving towards each other [2]. We have considered the exact 
solution known as the breather to serve as a benchmark to our numerical algorithm. This 
solution describes the interaction of a pair of kink/anti-kink interacting to form a periodic 
bound state, and is given by. 


(p{t, x) = 4 arctan 


Vi — a;2 sin{oj{t + to)) 

oj cosh(\/l — ui'^x) 


( 2 . 12 ) 


where U! is the frequency of the breather. We have determined the initial data (2.9) corre¬ 
sponding to the above solution and evolve the dynamical equations (2.6) and (2.7). Instead 
of comparing the exact solution with the approximate solution given by (2.12), we have 
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used the deviation from the conserved total energy E as the error control. More specifi¬ 
cally, after evaluating the total energy associated to the numerical solution at each instant, 
E{t), we determined the quantity, 

SE{t) = ^exact-^(t) ^ ^2.13) 

Eexa,ct 

where inexact is the conserved exact total energy. We have set oj = 0.5, to = 1-0 and evolved 
numerically the KG equation from t = 0 to t = 70 for increasing truncation orders, namely 
N = 10, 20, ..,50. In order to present a clear effect of increasing the truncation order, we 
have plotted the rms of 6E; 6E^^s, in terms of N. The result is presented in Fig. 2 with the 
expected exponential decay of <5Firms until the saturation due to round-off error is achieved 
for > 50. 

We have performed an additional numerical test of checking the conserved total energy 
of the collision of a pair of kink-antikink in the (/>^ model {p = 1). The initial configuration is 
given by uq) = —1 -|- tanh “ tanh > where we have set xq = 12 and 

fixed the boost parameter as u = 0.3. We have evolved the field equations (2.6) and (2.7) 
with the truncation orders N = 50,100,150, 200 and evaluated, in each case the error rms 
in the energy conservation from t = 0 to t = 100 when both kinks have collided and moved 
apart each other. Fig. 3 depicts the results. Accordingly, even for a modest truncation 
order of = 50, the energy is preserved to about a one part in 10^, whereas for N = 200 
the conservation of energy achieves about less than one part in 10^. 

3 Collision of two-kinks defects 



Figure 4. Illustration of the initial distribution of the energy density p{x, 0) corresponding to a 
pair of a two-kinks with impact velocity u = 0.2 described by Eq. (3.1). The two-kinks are placed 
at Xq = 15, and this initial distribution is similar to four kinks of (f)'^ model. 

As we have mentioned in the Introduction, the collision of a kink and anti-kink in 
the model exhibits a complex structure that depends on the impact velocity. This 
complex structure is expressed by the fractality associated to the transition between two 
possible outcomes: the formation of a bounded oscillating configuration and reflection of 
the kinks after the collision. These outcomes arise when the impact velocity is smaller or 
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Figure 5. Three-dimensional plot of the energy density for the collision of two-kinks defects for 
u = 0.3. Both two-kinks interacts in a very intricate way resulting in their complete disruption. 
The scalar field is radiated away leaving almost no structure about the origin 


greater than a critical velocity, Vc, respectively. However, for some values of the impact 
velocity below Vc the kinks become trapped, but it appears an intricate pattern of resonant 
windows in which they can escape to inhnity. The accepted explanation for this feature 
is a nonlinear competition or transfer of energy between the translational and vibrational 
modes of individual kinks [6]. 

We describe here the head-on collision of a two-kink and two-antikink pair for p > 3. 
The starting point is to establish the initial data representing both two-kinks configurations 
moving initially toward each other with velocity u. The initial data is given by, 

0) = -1 -k (pxix + xo, 0) - (I)k{x - xo, 0), (3.1) 

where we have chosen xq as the initial distance of the two-kinks defects. It is interesting to 
point out that the collision of both two-kinks might seem as the collision of four standard 
kinks. We show in Fig. 4 the initial profile of the energy density p{x,0). 

We begin with p = 3 that represents the most immediate generalization of the (j)‘^ model 
[p = 1), however most of the main features are also found on p > 3. We have performed 
numerical experiments with the velocity u as the free parameter. In all simulations, we have 
considered N = 250 and taking into account only the even modes in the spectral expansion 
(1.3) due to the symmetry of the problem. We have also monitored the error measured by 
the relative deviation of the total energy given by Eq. (2.13). It is worth mentioning that 
the error is sensitive to the choice of the map parameter Lq (cf. (2.2)). Thus, after some 
trial and error we have fixed Lq = 90 such that the error does not surpass 0.01%. We have 
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considered acceptable, in spite of being much greater than the typical 10“®% obtained for 
the collision in the (j)'^ model. We attribute this discrepancy to the nature of the outcomes 
arising after which we are going to describe. 



Figure 6. Collision of two-kinks defects for u = 0.5. Notice the formation of an oscillon at the 
origin and two symmetric oscillons that collide after bouncing. The outcome is an oscillon at the 
origin. 

The initial impact velocity u is the free parameter that plays a central role in the 
dynamics of the collision. We have illustrated in Fig. 5 the result of the collision for u = 0.3 
with the three-dimensional plot of the energy density p^j, given by Eq. (1.5). Both two-kinks 
defects collide in a very intricate way resulting in their complete disruption radiating away 
almost all scalar held. Eventually, a slight amplitude oscillatory bound structure about the 
origin survives. We have interpreted this bounded state as an oscillon-like structure. We 
remark that the observed fragmentation of both two-kinks in the hrst stages of the collision 
is a general feature regardless the value of the impact velocity and the value of p > 3. From 
the computational perspective, the dispersion of a large amount of the scalar held requires 
collocation points to considerable distances and, for this reason, large map parameter. We 
remark that this aspect is the leading cause of the relatively high error if compared with 
the collision in the model. In this last model, most of the scalar held is trapped at the 
origin or escape as traveling waves in localized distributions of energy. 

We have noticed that by setting higher impact velocities a complete fragmentation 
of both two-kinks. However, the fragmentation gives rise to the appearance of moving 
symmetric localized distributions of energy, together with a small and oscillating structure 
about the origin (cf. Fig. 5). We have understood these moving localized distributions 
energy as oscillon-like structures, from now on we call them oscillons. The formed oscillons 



Figure 7. Illustration of one of the critical configurations obtained for u = 0.1226. After the 
collision, two prominent symmetric oscillons remain approximately at rest after receding to a certain 
distance. We have noticed the formation of an oscillon at the origin. Eventually, both symmetric 
oscillons coalesce at the origin. 


recede to a certain distance before returning to collide to form a small amplitude oscillon 
about the origin. The effect of increasing the impact velocity is to increase the distance both 
oscillons move apart from each other before bouncing. In particular, these two symmetric 
prominent oscillons resemble the energy densities of the standard kinks of the (j)'^ model. In 
Fig. 6 we show the collision for u = 0.5, where the presence of these structures is present. 
After bouncing the oscillons collide resulting in a small amplitude remnant located at the 
origin survives. 

We have identified several configurations that appear between the previous outcomes, 
namely dispersion followed by an oscillon at the origin, and the formation of moving oscil¬ 
lons. These new structures emerge after fine-tuning the impact velocity to some particular 
values. We have named them critical configurations characterized when two symmetric 
oscillons remain at rest after moving apart each other to some distance. The time the 
oscillons are in this quasi-stationary phase increases if the impact velocity approaches one 
of its critical values. We have provided an illustration of one of the critical behaviors in 
Fig. 7 with u = 0.1226. It becomes evident the formation of three main oscillons after the 
collision: one at the origin and two symmetric. These last oscillons move to some distance 
and stay at rest for some time. Eventually, they coalesce at the origin. 

We can further fine-tune the impact velocity inside another interval such that the 
distance the oscillons recede becomes smaller. As before, the time both oscillons remain 
frozen depends upon which the fine-tuning to the critical velocity. We have presented in Fig. 
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8 two illustrations of the critical configuration with the three-dimensional plots of the scalar 
field and the energy density. In this case the impact velocity u = 0.39582. The interaction 
between both two-kinks occurs at f ~ 60, where two oscillons emerge, ft becomes evident 
that the oscillons remain at a fixed position forming a stationary structure, after moving 
apart each other. This phase lasts until t ~ 200. 



Figure 8. Critical configuration formed after the collision of the two-kinks with impact velocity 
u = 0.39582. In the upper plots, we have presented the three-dimensional plots of the scalar field 
and energy density. The stationary configuration is formed at t r; 60 and lasts to t « 200. In the 
lower plots, we have depicted the scalar field and energy density profiles evaluated at t = 180 (blue 
lines) together with the exact profiles given by Eqs. (3.2) and (3.3) (black circles). The agreement 
between the exact and numerical profiles is excellent. 

We are compelled to figure out the nature of the stationary configuration shown in 
Fig. 8. In this direction, we have considered the scalar field profile at t = 180 during its 
quasi-stationary phase. The objective is to find an exact solution that reproduces the scalar 
field and the energy density profiles. We have found that the best candidate that reproduce 
these numerical patterns is given by, 

(j)e{x) = (f)Q + Aq tanh'^ {box), (3-2) 

where </>Oj^o ^-^d bo are arbitrary constants and q is an even number. The corresponding 
energy density reads as. 
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(3.3) 
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- 1/9 




where 64>e = (f^e — 4>0- The best fit of the numerical scalar field profile at t = 180 as shown 
in Fig. 8 is obtained with the following parameters: (/>o = 0, Aq = —0.9833, 6o = 0.340 and 
q = 8. By inserting these parameters into the energy density (3.3) we have also reproduced 
the numerical profile of the energy density (cf. Fig. 8). 






Figure 9. Numerical profiles of the scalar field and energy density (blue lines) of the critical 
configuration together with the corresponding exact profiles given by Eqs. (3.2) and (3.3) (black 
circles). From up to down the numerical profiles result after the collision of a two-kink pair in the 
cases of p = 3 at t = 135.8, p = 5 at t = 350 and p = 7 at f = 400, respectively. The impact 
velocities are u = 0.1150, 0.11865 and u = 0.160. The lumplike critical solutions have the following 
parameters: (i) 0o —= —0.06656, Aq ~ —0.8666, bo ~ 0.3333 and g = 2; (ii) (j)o = 0, Aq ~ —0.9667, 
bo = 0.20 and q = 32; and (hi) (j)o ~ 4.8 x 10“^, Ao ~ —1.0867, bo — 0.1466 and q = 12. 

We remark that the solution (3.2) is a slight generalization of the static solution (1.3). 
As we have mentioned the case q = 2 represents an unstable topological lump-like defect. 
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However, one can verify that the remaining cases, namely q = 4,6,8..., possess the same 
properties of a lump-like defect [27, 28], and consequently we may consider these cases as 
representing lump-like defects. The numerical experiments indicate that a lump-like defect 
with q = 8 emerges as a critical configuration after the collision of two topological two-kinks 
defects [p = 3) with impact velocity u = 0.39582. As far as we are concerned, it is the 
hrst time an unstable topological lump-like defect is formed as the result of a collision of a 
topological defect pair. 

Further numerical experiments have shown the appearance of critical configurations 
in several intervals of the impact velocity. As a consequence, we conjecture that they are 
indeed lump-like solutions given by Eq. (3.2) with distinct values of q. For instance, by 
adjusting the impact velocity to u = 0.1150, we have obtained that the critical configuration 
is approximately represented by a lump-like solution with q = 2 (cf. Fig. 9). We have 
noticed in this case that the critical solution survives during a small interval (from t ~ 125 
to t ~ 150). Also, the value of 4>q is distinct from zero and can vary with time. It means 
that the critical configuration can oscillate about the exact static solution (3.2). 

We have found critical configurations when a two-kinks pair collide in the case p > 3. 
After settings = 5 (in this case xq = 20) and impact velocity u = 0.11865, we the lump-like 
solution has q = 32 and reproduces the critical configuration. Moreover, if p = 7 (xq = 30) 
we noticed that a critical configuration emerges when u = 0.160. We have presented the 
corresponding numerical and exact profiles of the scalar field and the energy density for 
these cases in Fig. 9. 

The appearance of the critical configurations can be the result of a nonlinear balance 
of energy transfer between the vibrational and translational modes of the two-kinks. The 
nonlinear equilibrium is achieved after an appropriate fine-tuning the impact velocity. In 
addition, the approach to the unstable lump-like solutions (3.2) indicates the existence 
of channels os attraction that is a new feature of these unstable solutions. Thus, these 
solutions might be saddle critical points in the abstract space of all possible solutions of 
the Klein-Gordon equation (1.1) with the potential given by Eq. (1.2). 

4 Discussion 

In this paper, we have studied the head-on collision of a new class of kinks known as two- 
kinks defects introduced by Bazeia et al [14]. We can understand the two-kinks as composed 
of two standard kinks of the cj)^ model. For this reason, the collision of a pair of two-kinks can 
be viewed as a collision of four standard kinks. The Klein-Gordon equation was integrated 
numerically with a simple and efficient algorithm based on the Gollocation method. We 
have presented numerical tests that confirm the accuracy and the rapid convergence of the 
algorithm. 

We have considered the collision of two-kinks defects for odd p > 1. In general, the 
outcome is the disruption of the two-kinks regardless the impact velocity. The scalar field 
is almost radiated away leaving behind a small amplitude oscillon at the origin. However, 
depending on the impact velocity there are two possibilities. The first is a direct formation 
of a small oscillon at the origin with the emission of almost all scalar field. The second is 
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the formation of outward moving oscillons that eventually bounce and collide to form the 
final oscillon at the origin. 

The new feature exhibited by the numerical experiments is the presence of several 
structures that interpolate the above intermediate outcomes. We have called these struc¬ 
tures as critical conhgurations. In these conhgurations, two symmetric oscillons remain at 
rest for some time after moving apart to each other. The time the oscillons evolve in this 
quasi-stationary phase depends on the hne-tuning of the critical velocity to some values. 
Most importantly, we have identified all critical configurations as belonging to a family of 
lump-like defects of topological nature and described by Eq. (3.2). It constitutes a new 
feature arising from the collision of a pair of topological defects, in this case, two-kinks 
defects with p > 3. 

We remark that a more thoroughly detailed numerical and analytical analysis are nec¬ 
essary for a complete understanding of the appearance of unstable lump-like defects. We 
have suggested that such a feature might be a consequence of a nonlinear equilibrium of 
the energy transfer between the vibrational and translational modes of the two-kinks. In 
this instance, the determination of the internal modes of two-kinks turns to be necessary. 
In Ref. [14] we found the zeroth mode, and the remaining are obtained after solving the 
following Schrodinger-like equation for the eigenmodes Xn{x), 

' (f ((fV\ 
dx"^ ^ V dctP' J K 

with an unusual property of the effective potential diverges at the origin. Another interest¬ 
ing consequence is the characterization of the family of solutions (3.2) as saddle-points in 
the abstract space of solutions of the Klein-Gordon equation (1.1). 
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